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Enteroviruses (EVs) are important human pathogens associated with various chnical syndromes. This study 
represents an overview of non-polio enteroviruses (NPEVs) isolated from acute flaccid paralysis (AFP) 
surveillance in Shandong Province, China from 1988 to 2013. Altogether 792 and 170 NPEV isolates were 
isolated from stool specimens of 9263 AFP cases and 1059 contacts, respectively. Complete VPl sequencing 
and typing on all 962 isolates revealed 53 NPEV types in which echovirus (E) 6 (7.6%), E14 (7.6%), Ell 
(7.4%), coxsackievirus (CV) B3 (7.4%), E25 (5.6%), CVB5 (4.9%), E7 (4.5%) and EV-A71 (4.4%) were the 
eight most commonly reported serotypes. Distinct summer-fall seasonality was observed, with June- 
October accounting for 79.3% of isolation from AFP cases with known month of specimen collection. 
Increase of isolation of EV-A71 and CVA — the predominant pathogens for the hand, foot, and mouth 
disease — was observed in recent years. Sequence analysis on VPl coding region of EV-A71 and E6 suggested 
Shandong strains had great genetic divergence with isolates from other countries. The results described in 
this study provide valuable information on the circulation and emergence of different EV types in the 
context of limited EV surveillance in China. 



Enteroviruses (EVs) {genus Enterovirus, family Picornaviridae) are common human pathogens. EV infection 
is usually asymptomatic, but sometimes it is associated with diverse clinical syndromes ranging from minor 
febrile illness to severe, potentially fatal diseases such as aseptic meningitis, encephalitis, paralysis, myo- 
carditis and neonatal enteroviral sepsis'"*. EV infections maintain to be an important public health problem. It is 
estimated that 10-15 million symptomatic enterovirus infections occurred in the United States each year''. 

EVs are small, nonenveloped, single-stranded RNA viruses. EV RNA is approximately 7.5 kb long and encodes 
a polyprotein that is processed to yield the mature structural (VPl to VP4) and nonstructural proteins'. EV 
serotypes have traditionally been classified into echoviruses (E), coxsackieviruses (CV) group A and B, and 
polioviruses (PV) by neutralization test, and molecular typing method based on VPl sequences has led to the 
discovery of a large number of new EV types'* ''. So far, human EVs comprises more than 100 types which are 
classified into 4 species, EV-A to EV-D'. 

PV infection is known to be associated with acute paralytic poliomyelitis. The standard approach recom- 
mended by World Health Organization (WHO) for polio surveillance is the detection and investigation of acute 
flaccid paralysis (AFP) cases, which includes standardized virological analysis of faecal samples of the patient, 
and/or sometimes those from contacts. Systematic AFP surveillance was initiated in China in 1994, and earlier 
experience on a small scale in Shandong Province could be dated back to 1988. Since 1996, annual reported AFP 
cases remained stable around 400. Besides PVs, a considerable proportion of non-polio enteroviruses (NPEVs) 
can be obtained from AFP surveillance as well. Since there is no specific EV surveillance system in mainland 
China yet, the AFP surveillance is the only data source for understanding EV circulations. 

In China, large-scale outbreaks of hand, foot, and mouth disease (HFMD) occurred repeatedly since 2007, and 
subgenogroup C4 was revealed to be the sole viral genetic lineage circulating in mainland China since then**. 
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Figure 1 | Annual numbers of AFP cases and contacts and isolation in Shandong Province, China, 1988-2013. 



However, existence of C2 subgenogroup of EV-A71 in Shandong in 
1996 was demonstrated in AFP surveillance', suggesting retrospect- 
ive investigation on EVs from AFP surveillance can provide valuable 
information on EV circulation. A study in Yunnan Province has 
identified 46 NPEV serotypes from AFP surveillance in 1997-2000 
and 2004'°. However, to the best of our knowledge, no similar study 
has been reported in other provinces of China yet. 

Shandong is a coastal province with a large population (—96 mil- 
lion) and major ports that could potentially serve as portals for 
importation of exogenous viruses. Here we describe the molecular 
typing and transmission patterns of all NPEV isolates detected dur- 
ing 26 years surveillance along with the phylogenetic analysis of two 
pre-dominant types EV-A71 and E6. 

Results 

Cases and NPEV isolation. The AFP surveillance was conducted in 
all 138 counties of Shandong Province and more than 600 sentinel 
hospitals were included. Hospitalized patients < 15 years of age were 



180 



reported, and male to female ratio was 1.9:1 for all cases from 1990 to 
2013. The data for the cases in AFP surveillance in the years 1988 and 
1989 were unavailable. Nevertheless, 4 NPEV strains were isolated 
during that period. During 1990-2013, a total of 9263 AFP cases were 
reported. EV isolation was performed on all stool specimens and 788 
NPEV strains were isolated from stool specimens of these cases. Also, 
specimens from 1059 contacts of AFP cases were collected during this 
period and 170 NPEVs were isolated. The annual numbers of reported 
AFP cases and isolates were illustrated in Fig. 1. AFP surveillance in 
China increased substantially throughout the 1990s, but the 
surveillance and report before 1995 was not so active, resulting in 
the low level of NPEV isolation. It should also be notified that the 
absence of information of contacts and related isolates before 2007 
(Fig. 1) are due to the unavailable information, not indicating that 
no specimens from contacts were collected at that time. 

The monthly distributions of AFP cases and NPEV isolation were 
illustrated in Fig. 2. AFP reports increased in summer and autumn 
months with a peak in July (1100 [11.9%] of9255 reports with known 
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Figure 2 | Numbers of AFP cases and NPEV isolates from AFP cases by month of specimen collection — Shandong, 1990-2013. 
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Table 1 | Annual isolation of different EV types from AFP cases in Shandong Province, China, 1 988-201 3 
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month during 1990-2013). However, a more prominent seasonality 
of NPEV detection from AFP cases was observed with June-October 
accounting for 79.3% (629 of 792) of isolation with known month of 
specimen collection. 

Serotype. Molecular typing was performed successfully on aU 962 
isolates from AFP cases and contacts during 1988-2013, and a total 
of 53 NPEV types were identified (Table 1, 2). Of these, the 16 most 



commonly reported serotypes (E6, E14, Ell, CVB3, E25, CVB5, E7, 
EV-A71, E3, E30, CVBl, CVA9, CVB2, CVA24, E13 and CVA4) 
accounted for 74.6% of all isolation. The eight most common 
serotypes accounted for approximately half of aU reports (49.3%). 
A total of 19 isolates belonging to 9 new EV types was identified 
including EV-A76, EV-A90, EV-B73, EV-B74, EV-B75, EV-B80, 
EV-B87, EV-B97 and EV-C96 (Table 1, 2). Co-infection of 
different serotypes was observed in 15 patients and 1 contact. 
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Table 2 | Annual isolation of different EV types from contacts of 
AFP cases in Shandong Province, China, 1 988-201 3 

Year 

Species Type 93 01 05 07 08 09 1 0 1 1 12 13 Total 
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During 1988-2013, NPEVs belonging to species A, B and C 
accounted for 10.5% (101 of 962), 84.4% (812 of 962) and 5.1% (49 
of 962) of total isolation, respectively. No EV-D virus was identified. 
EV-B is the major species in most years of the study period. However, 
dramatic increase in EV-A constitution in AFP cases was observed in 
recent years (Table 1). Especially in 201 1 and 2013, EV-A accounted 
for 39.0% (16 of 41) and 56.3% (18 of 32) of annual isolation, respect- 
ively. Common HFMD pathogens, such as EV-A71, CVA2, CVA4, 
CVA6, CVAIO and CVA16 were all present in the 26-y surveillance. 
Interestingly, several EV-A viruses had been observed in AFP sur- 
veillance before the onset of national HFMD epidemic in China in 
2007. The first isolation of EV-A71, CVA2, CVA4, CVA6 and 
CVAIO can be dated back to the year 1996, 1996, 1996, 1992 and 
2004, respectively. 

The most commonly isolated serotypes changed over time. Long- 
term trends of circulation of eight predominant NPEV serotypes 
were illustrated in Fig. 3, and epidemic pattern of circulation was 
observed for all serotypes including EV-A71, which remained rare 
before 2010 whereas increased dramatically since 2011. In 2011 and 
2013, EV-A71 accounted for 14.6% (6 of 41) and 31.2% (10 of 32) of 
annual isolation, respectively. Epidemic pattern of circulation of E6, 
El 1 and CVB5 was most distinct with regular sharp increases every 
3-6 years, which usually last for 1-3 years (Fig. 3). No geographical 
aggregation was observed for these NPEV types. 

VPl sequence analysis of EV-A71. Altogether 42 EV-A71 isolates, 
32 from AFP patients and the rest 10 from contacts, were identified 
during the AFP surveUIance in 1988-2013. During the period from 
1988 to 2002, no EV-A71 strains were isolated except for one C2 
subgenogroup strain in 1996'. Since 2003, continuous detection of 



EV-A71 was observed and they all belonged to C4a subgenogroup. 
Homologous comparison revealed 91.8% to 99.8% VPl nucleotide 
sequence similarities among C4 subgenogroup strains from AFP 
surveUIance, and 81.4% to 83.6% with prototype strain BrCr. Phylo- 
genetic analysis was performed on the complete VPl sequences of 
the isolates from AFP surveillance with global reference sequences 
and those previously obtained from patients with HFMD and aseptic 
meningitis in Shandong Province" (Fig. 4). Except for the C2 
subgenogroup strain in 1996, Shandong AFP strains were closely 
related to local HFMD and aseptic meningitis isolates from 2007 
to 2010 and belonged to C4a subgenogroup. 

VPl sequence analysis of E6. A total of 73 E6 strains were isolated 
during the AFP surveUIance in 1988-2013. Homologous comparison 
revealed 77.3% to 100% VPl nucleotide sequence similarities among 
themselves and the overall mean p-distance value is 0.162. Phylo- 
genetic analysis was conducted on the complete VPl sequences of the 
isolates with global reference E6 sequences. Shandong E6 strains 
segregated into 4 clusters, A-D (Fig. 5). Cluster A consisted of iso- 
lates from 2007 to 2012, indicating their circulation in Shandong is 
recently established. Nine strains from other provinces of China also 
belonged to this cluster. Cluster B consisted of isolates from 1988 to 
2002. Members of this cluster disappeared for 11 years, suggesting it 
might become extinct in Shandong. Cluster C was composed of 
isolates from 1993 to 2011, indicating this cluster had circulated in 
Shandong province for a long time. Great genetic divergence was 
observed with other E6 strains (mean p-distance between groups, 
0.204). Cluster D consisted of only 1 strain from 2007. It has close 
relationship with a strain from aseptic meningitis patient in Anhui 
Province, China in 2005. They were located in a branch consisted of 
strains from Australia, France, Russia, and Korea. However, no other 
Chinese strain of this branch was observed, suggesting this cluster 
might be imported into China recently and circulated at a low level. 

Discussion 

Since the Global Polio Eradication Initiative (GPEI) launched in 
1988, the global incidence of poliomyelitis has dropped a lot. No 
cases due to wild-type poliovirus 2 (WPV2) have been identified 
since 1999'^. However, AFP surveillance is stUl of great importance 
in monitoring re-emergence of WPV poliomyelitis in many prev- 
iously polio-free countries and emergence of vaccine-derived polio- 
viruses (VDPVs)'"*. During AFP surveillance in Shandong, no WPV 
has been isolated since 1992, and VPDV had been detected in 2007 
(PI, 9 nt;Pl, 13 nt), 2009 (P2, 11 nt) and 2011 (PI, 10 nt), respect- 
ively. They were all ambiguous VPDV (aVDPV) which were clinical 
isolates from persons with no known immunodeficiency. In USA, the 
National Enterovirus SurveUIance System (NESS) has provided valu- 
able information for investigating temporal patterns of circulation of 
different serotypes, guiding outbreak investigations, and identifying 
targets for development of diagnostic assays and antivirals^. Whereas 
in China, EV surveillance based on human specimens is very limited 
and AFP surveillance provided the only information on EV 
circulation. 

Since the correlation of serotype with VPl sequence was demon- 
strated^, molecular techniques of EV typing are becoming increas- 
ingly avaUable and new EV types continue to be identified'*"*. In the 
26-y surveillance, a total of 9 new EV types were identified. The 
numbers of isolates of each new type ranged from 1 to 4, accounting 
for a small proportion of total isolation (2.0%, 19/962). According to 
a previous 5-y study in Yunnan Province'", some of these new EV 
types such as EV-A76, EV-B75, EV-B80 and EV-C96 were also 
found. There is a long geographic distance (>1800 km) between 
the two provinces, and the co-detection of these serotypes suggest 
they had a wide spread in mainland China, although with low levels 
of activity. However, the eight most common serotypes in the 
Yunnan study were E13, E14, E12, CVB3, E2, E6, EV-C96 and El, 
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Figure 3 | Annual number ofisolationofS predominant serotypes from AFP cases, Shandong, 1990-2013. Annual isolation ofE6,Ell,E14andEV-A71 
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which are quite different with the results of this study (E6, E14, Ell, 
CVB3, E25, CVB5, E7 and EV-A71). Considering the results of Ell, 
E13, E7, E6, CVA24, ESQ, E19 and E29 as the eight most common 
serotypes in Philippines in a 17-y study'^, the variation in consti- 
tution of most common serotypes suggests EV serotypes have in- 
congruent epidemic magnitude under different spatio-temporal 
conditions. 

In mainland China, nationwide epidemic of HFMD has occurred 
repeatedly since 2007. EV-A71 and CVA16 are the two major cau- 
sative agents in different areas, and other coxsackieviruses of species 
A, such as CVA2, CVA4, CVA6, CVAIO, CVA12, etc, have also been 
demonstrated to be involved^""^^. Consistently, an increase of isola- 
tion of EV-A71 and other EV-A serotypes from AFP cases was also 
observed (Table 1, Fig. 3). EV-A71 has been demonstrated to be 
associated with neurological diseases in young children including 
aseptic meningitis^', poliomyelitis-like paralysis^'*'^', and fatal 
brain-stem encephalitis^''. The increase in isolation from AFP cases 
suggests since the large-scale epidemic of HFMD in 2007, EV-A71 
associated serious neurological diseases have become a major public 
health concern in mainland China, which is consistent with the 
observation of emerging EV-A71 detection in aseptic meningitis 
patients in a previous study in Shandong Province^'. 

Shandong Province is located in temperate zone. EV infections are 
peaked in summer-fall seasons in temperate climates, and the results 
from oiu' study are consistent with the seasonality (Fig. 2). According 
to results from EV surveillance in USA in 1970-2005, long-term 



trends of circulation of different EV serotypes have two major pat- 
terns: epidemic and endemic. Epidemic pattern is characterized by 
substantial fluctuations in circulation levels over time, including 
large peaks when the serotype was among the most prevalent enter- 
oviruses reported for a given year. Endemic pattern has stable and 
usually low-levels of circulation with few distinct peaks''. In this 
study, all eight most common serotypes displayed epidemic pattern 
of long-term trends, although these data are from AFP surveillance, 
not specific EV surveillance system. EV-A71 had large-scale nation- 
wide epidemic in current China, and an increase of isolation was 
observed in AFP cases in this study, whereas no similar trend was 
observed in the surveillance from USA^'^". Further surveillance will 
help understand the long-term trend of EV-A71 circulation in China 
under the circumstance of continuous epidemics of HFMD. In main- 
land China, C4 was demonstrated to be the predominant subge- 
nogroup responsible for most documented HFMD outbreaks and 
epidemic since 1998*, and Shandong strains from meningitis patients 
in 2010 also belonged to subgenogroup C4a". VPl sequence analysis 
in this study revealed close genetic relationship among EV-A71 C4 
subgenogroup strains from AFP, HFMD and aseptic meningitis 
patients, suggesting C4 is an important pathogen of multiple clinical 
manifestations. 

E6 is one of the most frequently isolated serotype in many regions 
of the world. Its infection includes aseptic meningitis, meningoence- 
phalitis, rashes, gastrointestinal illnesses, neonatal hepatitis and 
pneumonitis. In China, outbreaks and sporadic aseptic meningitis 
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Figure 5 | Phylogenetic tree based on E6 VPl sequences of Shandong isolates and global strains. Red triangles indicate Shandong strains from AFP 
surveillance, and green squares indicate strains from other provinces of China. Circle indicates the prototype strain and VPl sequence of E29 strain JV-10 
serves as an outgroup. 
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Table 3 | Primers used for PCR amplification of the VPl region of 
EV-C strains 



Primer 


Sequence (5' to 3') 


Gene 


Position* 


UFl 


GATGAYTWYACIGMIGGIGG 


VP3 


2309-2328 


UF2 


ATGKIYATIYTIGSinYGT 


VPS 


2387-2406 


URl 


CCRTCITARAARTGISWRTAIGC 


VPl 


31 1 1-3089 


UR2 


YTGCCAIGTRTARTCRTCCC 


VPl 


3007-2988 


040= 


ATGTAYRTICCIMCIGGIGC 


VPl 


2951-2970 


DRl 


GCYWTRnYTGRTGICCRAA 


2A 


3408-3389 



*Nucleotide sequence coordinates are given relative to the genome sequence of PVl-Sabin 
(AYl 842 1 9) for orientation only. 



cases caused by E6 have been reported recendy^^'^"^. In the VPl phylo- 
genetic tree, Shandong strains were grouped into four clusters. 
Except that the only member of cluster D had close relationship with 
strains from France, Australia, Russia and South Korea, the rest 3 
clusters contained Chinese strains exclusively, suggesting the exist- 
ence of geographic barrier between mainland China and other 
regions of the world. Frequent travel might increase importation 
or exportation of different E6 lineages. 

Co-infection was detected in 15 patients and 1 contact. It was 
identified in two kinds of situations. First, in this study, entire VPl 
sequence of an EV strain was obtained via combination of 5 ' end and 
3' end sequences. So, if molecular typing on 5' end and 3' end PCR 
sequences revealed different types, this indicates a mixture of differ- 
ent EV types in one isolate. Second, RD and HEp-2 cell lines have 
been demonstrated to possess different sensitivity to different sero- 
types. RD is sensitive to echovirus and EV-A CVAs, while HEp-2 is 
sensitive to CVBs and some EV-C types. So, isolation of different EV 
types from a single specimen was observed from different cell lines. 
In comprehending the surveillance results of this study, the limita- 
tion of virus isolation should be considered. Some serotypes, such as 
some coxsackie A viruses of EV-C species, cannot produce visible 
CPE in these cell lines. So, these viruses cannot be recovered in the 
cell culture method. Hence, the actual EV distribution in the popu- 
lation might not be completely identical with that reflected by this 
study. Also, this study describes EV typing results from AFP surveil- 
lance. Since some NPEVs, such as EV-A71, had been demonstrated 
to be the pathogens of AFP, the results of this study might be different 
from usual EV surveillance (infections often milder) and do not 
necessarily directly reflect to general EV circulation. 

China has a vast territory and a large population. The EV circula- 
tion in Shandong Province is not necessarily consistent with that of 
other regions in China. However, EV has the ability of transmission 
over long distance, and frequent transmissions over cities or pro- 
vinces in China have been observed previously**'^''. So, the results of 
serotype spectrum, long-term pattern of circulation and molecular 
epidemiology can serve as a reflection of the situation in other parts 
of Eastern China. Further similar studies from other regions wfll help 
understand the EV circulation in the country. 

Methods 

Ethics statement. Ethical approval was given by the Ethics Review Committee of 
Shandong Center for Disease Control and Prevention, and the methods were carried 
out in accordance with the principles of the Declaration of Helsinki. Written 
informed consents for the use of their clinical samples were obtained from all subjects 
(the legal guardians of the patients and contacts). 

Virus isolation. Stool samples from the AFP patients and contacts were collected and 
processed according to standard procedures recommended by the WHO^\ Contacts 
were collected for AFP cases with (1) <5 years of age, (2) <3 doses of OPV 
immunization or unknown OPV history, (3) unqualified stool specimen, or (4) 
clinical diagnosis of poliomyelitis. RD (human rhabdomyosarcoma cell line), HEp-2 
(human epidermoid carcinoma cell line) and L20B (mouse cell line expressing the 
gene for the human cellular receptor for poliovirus) cell lines were used for virus 
isolation^\ All cell lines were gifts from the WHO Global Poliovirus Specialized 
Laboratory in the USA and were all originally purchased from the American Type 



Culture Collection. A total of 200 |.il of chloroform- treated stool solution was added 
to each of the cell culture tubes, and the inoculated cells were examined daily. After 
7 d, the tubes were frozen, thawed, and re-passaged, and another 7-d examination 
was performed. Infected cell cultures were harvested and used for further 
examination untU complete cytopathic effect was obtained. To avoid cross 
contamination, cell tubes of normal L20B, RD and HEp-2 cells served as negative 
controls. Isolates from RD or HEp-2 cell lines were re-passaged to L20B cell line, and 
were designated as NPEV if no CPE was observed. 

VPl amplitication, sequencing and molecular typing. Total RNA was extracted 
from 140 |,d of the infected cell culture using QIAamp viral RNA mini kit (Qiagen, 
Valencia, CA, USA) and reverse transcription- PCR (RT-PCR) was performed using 
Access RT-PCR System (Promega, USA) according to the manufacturers' procedures. 

Primer pairs 486/488, 487/489 and 040/01 1 for EV-A strains, and primer pairs 008/ 
013, 490/492, 187/011, 491/493 and 012/011 for EV-B strains were used for amp- 
lification of partial VPl sequences*-^-^^. As for EV-C VPl amplification, previously 
designed primer pairs 494/496, 495/497 and 040/011'*'^^ are not always efficient. So, 
we designed primer pairs UFl/URl, UF2/UR2 and 040/DRl (Table 3) in this study 
based on alignment of all available EV-C complete genome sequences from GenBank. 
These primers turned out to be efficient not only for all EV-C strains described in this 
study (CVA17, CVA20, CVA21, CVA24 and EV-C96), but also for Sabin PVl-3 and 
other available CVA13 and EV-C99 strains (from healthy population) of EV-C in our 
lab. Nevertheless, due to limited EVs in our lab, their performance on other EV-C 
viruses beyond the listed 10 serotypes had not been testified. Combination of 5' end 
and 3' end sequences yielded entire VPl coding sequences. As the serotype and 
species were unknown before typing, EV-B primers were first used for RT-PCR 
amplification. And EV-A and EV-C primers were used in case of negative results by 
using EV-B primers. In order to avoid cross contamination, a RT-PCR reaction using 
the RNA extracted from normal RD cell served as a blank control, and a negative 
control containing all the components of the reaction except for the template was also 
included. 

PCR positive products were purified and sequenced bi-directionally with the 
BigDye Terminator v3.0 Cycle Sequencing kit (Applied Biosystems, Foster City, CA), 
and sequences were analyzed by ABI 3130 genetic analyzer (Applied Biosystems, 
Hitachi, Japan). Molecular typing based on VPl sequences was performed using 
online Enterovirus Genotyping Tool version O.P^. 

Homologous comparison and phylogenetic analysis. Sequence analysis was 
performed on two predominant serotypes, EV-A71 and E6. The entire list of VPl 
sequences of global reference E6 and EV-A71 strains used in sequence analysis can be 
found as Supplementary Table Si online. Nucleotide sequence alignments were 
carried out by BioEdit 7.0.5.3 software'*. Phylogenetic trees were constructed by Mega 
4.0 using neighbor- joining method after estimation of genetic distance using the 
Kimura two-parameter method'^. A bootstrapping test was performed with 1,000 
duplicates. 
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